Structural stability and magnetism of FeN from first principles 
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Abstract 

In the framework of density functional theory (DFT), the structural and magnetic properties 
of FeN mono nitride have been investigated using the all electrons augmented spherical wave 
method (ASW) with a generalized gradient GGA functional for treating the effects of exchange 
and correlation. Calculation of the energy versus volume in hypothetic rocksalt (RS), zinc-blende 
(ZB) and wurtzite (W) types structures shows that the RS-type structure is more stable than the 
others. Spin polarized calculation results at equilibrium volume indicate that the ground state of 
RS-FeN is ferromagnetic with a high moment, while ZB-FeN and W-FeN are non magnetic. The 
influence of distortions on the stability is taken into account by considering FeN in two different 
face-centred tetragonal structures (fct): fct-rocksalt and fct-zincblende. The magnetovolume effects 
with respect to Slater-Pauling-Friedel model are discussed. The electronic structures analyzed from 
site and spin projected density of states are reported. A discussion of the structural and magnetic 
properties of FeN is given with respect to N local environment of Fe. 

PACS numbers: 71.15.Mb, 71.15.Nc, 71.20-b, 75.10.Lp, 74.25.Ha, 73.43. Cd 
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I. INTRODUCTION 



The study of the iron-nitrogen system has attracted much scientific importance for basic 
research as well as for technology. It was extensively investigated over the two decades 
(1950-1970) following its first discovery by Jack [l|. A large number of results have been 
reported for Fe-rich nitrides, see ^ and references therein. This was motivated by their 
potential appUcatio,. such as pig„.e„ts fo. high density .nagnet.c tecotd™ g On the 
contrary, only few reports are available in the case of the mononitride FeN [J, |5| , which is 
considered as interesting material in the emergent spintronic field (sources of spin injection 
to semiconductors or diluted magnetic semiconductors). Recently, thin films of FeN have 
been synthesized and two possible crystal structures were reported One is 

the rocksalt (RS) structure with a lattice constant of a = 8.638 bohr, and by performing 
^^Fe Mossbauer spectroscopy measurements the authors suggested rocksalt FeN to be an 



ant if erromagnet 



,|7|]. The other structure is of zinc-blende (ZBj-tjpe with a lattice constant 
of a = 8.185 bohr and a micromagnetic character 8|]. 

In order to achieve a better understanding of their physical properties and to discuss the 
experimental discrepancies, some theoretical investigations have been undertaken. From the 
available literature a clear controversy exists between the studies carried out in the frame- 
work of density functional theory. Shimizu et al. [l^ have concluded that a ferromagnetic 
rocksalt FeN is more stable than non magnetic zinc-blende nitride. The same results were 
reported later on the stability of the^FeN RS-type structure, but without a complete agree- 

12 1 . However, Lukashev et al. [l3| have joined some 



ment for the magnetic ordering |ll, 

others [l^ in predicting a stable non magnetic ZB structure for FeN. In this context and 
in order to provide an improved theoretical explanation, we carried out new comprehensive 
spin density functional calculations in both ZB and RS phases. Further, for a better in- 
vestigation of the magnetic properties with respect to the structural aspects, we have also 
considered FeN in the hexagonal wurtzite-type structure. With the sophisticated techniques 
of synthesis reached recently (synthesis of thin films and nanowires, ...), the structural 
distortions may have important effects on the investigated physical properties. To take into 
account the effects of some possible distortions, we have also investigated FeN in the face 
centered tetragonal fct structure. As a matter of fact, the manganese mononitride (MnN) 
has been recently found in (fct-RS) -type structure i.e., a tetragonally distorted rocksalt 
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structure [15|]. In our case, we have considered FeN in both tetragonally distorted rocksalt 
structure (fct-RS) and tetragonally distorted zinc-blende (fct-ZB) structure for sake of com- 
pleteness. To the best of our knowledge, this is the first theoretical study of FeN in wurtzite, 
fct-RS and fct-ZB type structures. The structural, electronic and magnetic properties are 
thus reviewed by analyzing the calculated electronic structure, total energies and crystal 
field influence. A description of the computational details is given in section 2. Our results 
for the calculated total energy and magnetic moment are presented in section 3, with a 
discussion of the magnetovolume effects. Section 4 is devoted to the Slater-Pauling-Friedel 
model as applied to iron nitrides. In section 5 the non magnetic results are discussed in 
the framework of the Stoner theory of band ferromagnetism. Section 6 details the electronic 
structure via partial density of states analysis. Finally a conclusion and a prospective for 
future developments are given in the last section. 



II. COMPUTATIONAL DETAILS 



Our first principles calculations are performed in the framework of the density functional 
theory using the generalized gradient approximation (GGA) with the Perdew, 

Burke and Ernzerhof PBE parameterization for the exchange and correlation [l^. This 
parameterization scheme was preferred over the LDA (local density approximation) in view 



Id, 



121 ]. Self-consistent calculations were 



of the results on iron nitrides cited in this work 
carried out using the scalar-relativistic augmented spherical wave method (ASW) [19|, l2( 
based on the atomic sphere approximation (ASA). In this method, the value of the GGA- 
PBE kappa was fixed to 0.8. In the ASW method, the wave function is expanded in atom- 
centered augmented spherical waves, which are Hankel functions and numerical solutions of 
Schrodinger's equation, respectively outside and inside the so-called augmentation spheres. 
In order to optimize the basis set, additional augmented spherical waves were placed at 
carefully selected interstitial sites. The choice of these sites as well as the augmentation 
radii were automatically determined using the sphere-geometry optimization algorithm 21|. 
The basis set consisted of Fe(4s, 4p, 3d) and N(2s, 2p) valence states. The Brillouin zone 
integrations were performed with an increasing number of k-points (16x16x16) in order to 
ensure convergence of the results with respect to the k-space grid. The convergence criterion 
is fixed to 0.001 mRy in the self-consistent procedure and charge difference A Q=10~^ 
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between two successive iterations. In order to establish a reference for the spin-polarized 
calculations we started with a set of spin-degenerate calculations. We note that this does 
not represent a paramagnetic situation which would require heavy computations involving 
large supercells with random spin orientations. It allows assigning a role to the orbitals 
responsible of the magnetic instability toward s pin polarization in a mean field analysis 



using the Stoner theory of band ferromagnetism 22|. 



III. RESULTS FOR THE TOTAL ENERGY AND THE MAGNETIC MOMENTS 

In a collective electrons scheme, such as the one used here, the magnetization arises from 
interband spin polarization, i.e., it is mediated by the electron gas. This is opposed to 
the localized electron moments where magnetization arises from intraband spin polariza- 
tion such as in oxide systems, especially insulating ones. From this, it is expected that the 
magnetovolume effects will be large in intermetallic and insertion alloys systems. The theo- 
retical equilibrium volume and the structural preference of the FeN compound are obtained 
by calculating the variation of the total energy versus the volume of the cubic cell in the 
rocksalt and zinc-blende structures. For wurtzite, /ci-rocksalt and /ci-zinc-blende, the same 
procedure (total energy versus the volume of the unit cell) was done twice, in order to obtain 
the two equilibrium lattice constants a and c which determine the equilibrium volume. The 
stability toward magnetism in each phase is given by comparing the spin-polarized (SP) and 
the non spin-polarized (NSP) total energy values at theoretical equilibrum volume. With 
the same k-points grid of the reciprocal lattice and for different pairs (E,V), after a fitting 



with a Birch equation of state [25|], we obtain after convergence a quadratic curve with a 
minimum (Eo,Vo). Figures 1-a and 1-b show the energy versus volume curves of FeN in 
the RS and ZB structures for NSP and SP cases, respectively. For the wurtzite structure, 
figure 1-c shows the energy variation versus volume at the equilibrium c lattice constant, i.e. 
the variation of energy versus the a lattice constant. Because the /ci-rocksalt and fct-zinc- 
blende structures are derived from the cubic ones, the same variation is observed and the 
results are not plotted here. The results for theoretical lattice constants, total energy and 
other equilibrium properties are summarized in table I for all the studied structures with 
the experimental data for the lattice constants. 

Figure 1-a shows that the ferromagnetic state of FeN in the RS phase is preferred at the cal- 
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culated equilibrium volume to the non magnetic state. On the contrary, the magnetic order 
in ZB and W phases can be favored but only at higher volume, see figure 1-b and figure 1-c. 
At equilibrium, the energy difference between the ferromagnetic and non magnetic states in 
these two structures are nearly zero, thus instability toward magnetism could exist. This 
instability may be the reason of the micromagnetic character observed experimentally for 
FeN in ZB structure. From table I our estimated theoretical lattice constant in the SP case 
(a = 8.219 bohr) agrees well with the experimental data for the ZB phase. However, a large 
difference exists between the obtained theoretical lattice constant (a = 7.769 bohr) and the 
experimental reports (ag^p = 8.638 bohr) in the RS structure. It is important to mention 
that our results for this point are in good agreement with all the theoretical studies of the 
FeN reported in the literature (even those predecting a ZB stable structure for FeN), but 
surprisingly, up to now, no one was able to reproduce the experimental data. 
The common explanation put forward to explain this large difference is that it would origi- 
nate firstly from surface effects because the samples were grown as thin films, and secondly 
from the nonstoichiometry of the elaborated samples. The values of the total energies in 
table I indicate that the FeN compound in the ferromagnetic state should be more stable in 
RS structure than in ZB or W ones. At equilibrium FeN in these two latter structures has 
a zero magnetic moment, while the unit cell of i^S-FeN possesses a moment of 2.62//^ (a 
moment of 2.51j»B per Fe-atom). Our results are in a good agreement with those reported 
by Shimizu et al. 10|] concerning the ferromagnetic order of FeN. A value of 2.67//B was 
reported Y.Kong but in a stable antiferromagnetic (AFM) i^S"- FeN 12|. For W-FeN, we 
think that this is the first study of this hypothetic structure, and thus a comparison of our 
results with others' cannot be made. 

We also performed antiferromagnetic (AF) calculations at the ferromagnetic equilibrium 
lattice constant in order to check for a possible instability toward such a magnetic con- 
:iguration for the ground state as suggested by Nagakawa er al. j^l and Hinomura et al. 
71]. Among possible spin arrangements we examinded for this purpose the [001]-alignement 
(AF-[001]) in which the spins are parallel within a layer but alternating direction between 
adjacents layers in [001] direction. From carefully converged calculations the results in Table 
I indicate that the energies of AF and ferromagnetic (FM) states for FeN in RS structure re- 
main very close. Moreover, in /ci-rocksalt structure, the former one seems more stable. The 
energy differences between AF and FM states are too small to make a clear conculsion on 
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the nature of the magnetic order in the two crystal varieties. On the other hand, the energy 
difference between RS-FeN and /ct-RS-FeN is very small, thus one can believe that under 
some elaboration conditions, FeN nitride can crystallize in a face-centred tetragonal-rocksalt 
structure. Because of this, the AF character observed experimentally can be explained either 
by this AF-FM instability in the i?S'-structure or by some lattice distortions which stabilize 

fci-rocksalt structure. 



the FeN compound in A 
In earlier investigations 26| on the magnetovolume effects in the Fe4N nitride, it has been 
shown that there is a transition from "low moment-low volume" to "high moment-large vol- 
ume" which resembles the moment versus volume dependence of 7Fe. This behavior seems 
to occur in all structures studied here of equiatomic FeN nitride. As it can be seen in figure 
2, the variation of the moment in i^S"- FeN (and in fct- RS-FeN which is not shown here) is 
very similar to those of Fell (located at the faces of the cubic cell) in Fe4N. In ZB-FeN (and 
in fct-ZB-FeN also not shown here) a sudden increase from zero to a value as high as 2.85//s 
occurs at high volume, like the behavior of Fel atoms (located in the corners) in Fe4N. The 
same variation of the magnetic moment is indicated for VK-FeN. 

An analysis of the effects of the Fe-N spacing on the magnetic moment value shows that 
the Fe-N distances, corresponding to equilibrium of zinc-blende, /ct-zinc-blende and wurtzite 
structures (see table I), give a non zero moment (around I/xb) in rocksalt or/and /ci-rocksalt 
structure. If we consider the similar increase of the magnetic moment value reported else- 



where in CsCl-type FeN 12|, one should recognize that at different atomic environements 
(crystalline structures) even with the same Fe-N distance, different values of the moment are 
obtained. Thus the major effect which prevails in the value of the magnetic moment is the 
nature of the d-p hybridization in the Fe-N bond which is related to the atomic environment 
(the number of the nearst-neighbors), and the distance Fe-N has only an average effect. 

IV. ANALYSIS WITHIN THE SLATER-PAULING-FRIEDEL MODEL 

The magnetic moment variation with valence electron count can be discussed by means 



of the Slater-Pauling (SP) model [2^. This is mainly illustrated by the SP curve which is a 
plot of the average magnetic moment /Jav with the valence electron count Zy for intermetallic 
alloy systems characterized by strong ferromagnetism. From the plot the variation of the 
average magnetization of an alloy with the solute concentration can be obtained. For most 
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of the magnetic systems, either the majority (|) or minority (j) spin population is known. 
The magnetization being provided by the difference of electron occupation between the 
majority (|) and the minoriry (|) spins, the following relationships can be established: m 
= nf - nj, and Z.^, = n| + n|. Hence m, the magnetic moment, can be obtained either 
as: m = 2 n{^) — Z^,, if n(|) is known, or m = Z^, — 2n(|), if n(|) is known. These 
two expressions for m describe two branches with opposite 45° slopes of the Slater-Pauling 
curve around which the experimental points are gathered. If one assumes that the magnetic 
moment mostly arises from d-band polarization, the quantity "magnetic valence" can be 
defined -here for systems where n^ f is known as Zm = 2nd^ — Z„. Upon alloying, the 



electron count of the d states changes in a discrete way as treated by Friedel [2J], so that 
rif^i is either (early transition elements) or 5 (late transition elements). For instance, for 
Fe, Z^, is 8 and Zm= 10-8 = 2, whereas for N (n^i = 0), = 3 and Zm = — 3 = —3. This 
leads to an alternative representation called Slater-Pauling-Friedel (SPF) curve where \iav 
is plotted against Z^,'- fJ^av = Zm + 0.6. In this writing , the figure 0.6 is a nearly constant 
contribution to the magnetic moment arising from s and p electrons. The calculated atom- 
averaged magnetic moments (fiav) via the expression above and from ASW calculations 
are summarized in Table II. The results of the present study (for FeN nitride) are shown 
in the last line of the table; the other data are from in which the SPF model was 



271 ] have 



applied for the whole series of the Fe-N systems. In their study, Matar and Mohn 
found a good agreement between SPF model and ASW calculation results, and all average 
moments decrease with increasing amount of nitrogen. While our obtained moments in 
ZB-FeN, fct-ZB-FeN and VK-FeN agree well with SPF model in finding a vanishingly small 
magnetization, a large discrepancy exists in the case of RS-FeN and fct-RS-FeN which carry 
large magnetizations. This interestingly points to the limits of the SPF model which is an 
average scheme dealing with the electronic configuration of the chemical species in presence, 
and does not account for the nature of the crystalline structure and the chemical bonding. 

V. MEAN FIELD ANALYSIS OF NON MAGNETIC CONFIGURATIONS 



Within the Stoner theory of band ferromagnetism 23|] which is a mean field approach, 
the large DOS at the Fermi level n(i?i?) is related to the instability of the non magnetic 
state with respect to the onset of intraband spin polarization. When n{EF)Is > 1 where 
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Is is the Stoner integral, calculated and tabulated by Janak [28|] for the elemental systems, 
then the non magnetic state is unstable with respect to the onset of magnetization. In the 
RS-FeN type with a density of states at the Fermi level n{Ep) = 71.596 states/Ry/Fe and 
a Stoner parameter Ig = 0.034 Ry the Stoner criterion is fulfilled (2.434). On the contrary, 
for the Z5-FeN and W-FeN we have niEp) = 14.106 states/Ry/Fe and niEp) = 15.152 
states/Ry/Fe respectively, the Stoner criterion (resp. 0.48 and 0.515) is not fulfilled. As 
a consequence, the magnetic behavior of FeN is consistent with the Stoner theory. As it 
was discussed in the literature, the magnetic moment of Fe atom in FeN [lO, Q, [l^ and 



generally in the Fe nitrides [2] is very sensitive to nearest-neighboring Fe-N distance and 
their d-p hybridization. Our results confirm that this latter is the most important effect. 



VI. DENSITY OF STATES IN THE FERROMAGNETIC CONFIGURATION 

Considering the results obtained for the electronic structure of FeN at equilibrium we dis- 
cuss the electronic structure from the plots of the density of states (DOS). Spin polarization 
results from a nearly rigid band shift between the majority spins (|) to lower energy and 
minority spin (i) to higher energy due to the gain of energy from exchange. The site and 
spin projected partial DOS (PDOS) are shown in fig. 3 for RS-FeN and ZB-FeN at the the- 
oretical equilibrium lattice constants. The VF-FeN PDOS which resemble those of Z5-FeN 
are not presented here. Energy low lying N(2s) states at ~ -10 eV as well as ES(PDOS) of 
vanishingly low intensity are not shown. As a matter of fact ES receive charge residues from 
actual atomic spheres and they are known to ensure for covalency effects in such systems 
Note that in RS-FeN the octahedral environment splits the Fe-3d states into t2g and eg man- 
ifolds, whereas in the tetrahedral environment of ZB-FeN, the d states are split into e and 
^2 manifolds. These respective environments result in totally different DOS features within 
the two structures: whereas they are continuous over the VB in i^S'-FeN more localisation is 
observed for the ZB-FeN. One consequence is that i?S'-FeN is a metal through the the itin- 
erant Cg states crossing the Fermi level while in ZB-FeN the system is close to an opening 
of gap at Ep. Within the valence band (VB) two parts can be considered. While the lower 
energy part, centered at -5 eV, is mainly composed of N2p states mixing with itinerant Fe-e^ 
ones, the higher energy part from -3 eV to 1 eV is dominated by the Fe-t2g states which are 
localized and thus responsible of the magnetic moment (RS-FeN). These states are moved 
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to lower energy for | spin population (upper panel) and hence are found below the Fermi 
level Ep, whereas for the minority-spin (j) the DOS peaks are located above Ep (lower 
DOS half panel). This is due to the exchange-splitting whose absence in the ZB-FeN PDOS 
causes mirror PDOS in upper and lower half panel of fig. 3b. Clearly the crystal field is 
unfavorable to the onset of magnetization in ZB-FeN and a fortiori in VF-FeN. 



VII. CONCLUSION 



With the use of the self-consistent DFT-based ASW method, we have investigated the 
structural, electronic and magnetic properties of the iron mononitride (FeN) in the rocksalt 
(RS), zinc-blende (ZB) and wurtzite (W) structures. We have also investigated the effects 
of distortion by considering FeN in a tetragonally distorted rocksalt fct-RS and in a tetrag- 
onally distorted zinc-blende fct-ZB structures. The calculated total-energy shows that FeN 
is stabilized in a ferromagnetic (FM)-rocksalt structure with a theoretical lattice constant a 
= 7.769 bohr and a magnetic moment of 2.62 ^ib- While ZB-FeN and W-FeN type seem to 
prefer a paramagnetic order at theoretical equilibrium, they possess nevertheless a nonzero 
magnetic moment at higher lattice constant. This involves strong magnetovolume effects 
which are equally accompanied by crystal field ones depending on whether the local environ- 
ment is octahedral or tetrahedral. The equilibrium energies of non magnetic and FM states 
of ZB-FeN are very close and this instability toward magnetism may be the reason of the 
micromagnetic character observed experimentally. We have also found a very small energy 
difference between RS-FeN and fct-RS-FeN, thus one can predict that in some special syn- 
thesis conditions such as the presence of strains in the growth of layers of iron mononitride, 
FeN could crystallize in the latter structure. The energies of AF and FM states in RS-FeN 
and fct-RS-FeN are too close, so it is difficult to make a clear conclusion on the nature of 
the magnetic order of the ground state for these structures. As it was suggested in the 
literature, the present work confirms that the magnetic properties of FeN are dominated by 
p — d hybridization of the Fe-N interaction. However, we have found that the Fe-N distance 
has only an average effect on the value of the magnetic moment. On the other hand, our 
results of FeN compound complete the finding on the competition between chemical and 



magnetovolume effects in the series of iron nitrides: FegN, Fe4N, FesN as well as Fe2N 



29|. 



For FeN in ZB and W-type structure a good agreement with Slater-Pauling-Friedel model. 
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connecting average magnetization with valence, is found; but when RS-FeN is considered a 
large difference exists so that other effects related with the nitrogen environment of Fe {Oh 
versus or C4„) should become prevailing. As a prospective we think of completing the 
explanation of the FeN properties in two complementary directions: 

• by considering slab and surface calculations, and 

• by taking into account non (i.e., sub- or super-) stoichiometries in extended lattices. 
This should comply with the nature of the samples elaborated experimentally which 
point to N substoichiometric FeNl x compositions. 

Such investigations will call for heavy calculations involving other computational frameworks 



30|. The studies 



such as the use of pseudopotential based codes for geometry optimizations 
are underway. 
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TABLE I: Calculated equilibrium properties (in atomic units(bohrs) lRy=13.6eV) for FeN nitride 
in five different structures: Rocksalt (RS) , zinc-blende (ZB), wurtzite (W), face centred tetragonal- 
rocksalt (fct-RS) and face centred tetragonal-zinc-blende f/c^-ZiJj 



equilibrium 


Lattice 


total 


dpe—N{bohr) 


nifnt (ur) 


-L VJ 1^ V-> A IJJ.V->1J 


constants ibobr i 








NSP-RS 


a=7.610 


-2655.026 


3.788 




SF-RS 


a=7.769 

a=7.559"^ 

a=7.939'' 


-2655.047 


3.788 


2.60 


AF-RS 




-2655.046 


3.890 


2.36 


NSF-ZB 


a=8 200 


-2654 997 


3.534 




SF -ZB 


a=8.219 
a=8.24" 
a=7.933'^ 
Exp: a=4.332'' 


-2654.998 


3.551 


0.00 


mF-w 


a=5.778 
c=9.637 


-2655.018 


3.606 




SF-W 


a=5 809 
c=9.637 


-2655 019 


3.608 


0.00 


NSF -fct-RS 


a=7.778 
c=7.580 


-2655.022 


3.869 




SF -fct-RS 


a=7.844 
c=7.655 


-2655.044 


3.924 


2.63 


AF -fct-RS 




-2655.045 


3.924 


2.48 


NSF -fct-ZB 


a=8.090 
c=8.015 


-2654.992 


3.523 




SF -fct-ZB 


a=8.128 
c=8.015 


-2654.993 


6.663 


0.00 
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TABLE II: Application of the Slater Pauming Friedel model of the magnetic valence to iron nitrides 
(see text). 



Iron nitride 






msppifJ-B) 

= Zra+0.6 




reglN 


1 A A 

/ .44 


1 A A 

i.44 


Z.U4 


Z.Ll 


Fe4N'^ 


7.0 


1.0 


1.60 


1.67 


FesN*^ 


6.75 


0.75 


1.35 


1.44 


Fe2N'^ 


6.33 
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FIG. 1: Spin polarized (SP) and non-spin polarized ( NSP) total energy versus volume in: a) 
RS-FeN, b) ZB-FeN and c) W-FeN, at equilibrium cq lattice constant for the latter. 
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FIG. 2: Magnetic moment versus volume in: a) present work: FeN in RS, ZB and W structures. 
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FIG. 3: Spin polarized total and partial density of states (PDOS) of a) i^S-FeN and b) ZB-FeH at 
equilibrium lattice constants in the ferromagnetic state. 
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